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& Abstract 
>> 

Steady-state turbulence is generated in a tank of water and the trajectories of particles forming a compressible system on the surface 
l Pn are tracked in time. The initial uniformly distributed floating particles coagulate and form a fractal distribution, a rare manifestation 
of a fractal object observable in real-space. The surface pattern reaches a steady state in approximately 1 s. Measurements are 
made of the fractal dimensions D q (t) (q = 1 to 6) of the floating particles starting with the uniform distribution D q (0) = 2 for 
Taylor Microscale Reynolds number Re\ ~ 160. Focus is on the the time-evolution of the correlation dimension D-2.(t) as the steady 
CO state is approached. This steady state is reached in several large eddy turnover times and does so at an exponential rate. 
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■H 1, Introduction 



When studying a dynamical system in phase space, one 
of several scenarios may occur. The simplest cases involve 
either a fixed point or a limit cycle (periodic orbit). The 
most interesting scenario unfolds when studying dissipative 
systems undergoing chaotic evolution. When a system is 
driven out of thermal equilibrium, a phase point which pre- 
viously would have visited all regions of phase space with 
almost equal probability, now spends most of its time in a 
limited region. There it develops a complex fractal topol- 
ogy with a non-trivial fractal dimension [13]. Fractal di- 
mensions have been measured extensively for simple maps 
(eg. the Henon map) [1,6] as well as for real world systems 
[2] . Typically, the fractal behavior of such systems is stud- 
ied in an asymptotic limit which ignores the transient be- 
havior [1]. In this paper, we show that measuring a fractal 
dimension in the transient state is an effective way to char- 
acterize the evolving topology of this particular system. 



In this experiment we study the transient evolution of a 
fractal topology in a laboratory setting, namely, the clus- 
tering of floating particles at the surface of a turbulent fluid. 
If neutrally bouyant, non-inertial particles are introduced 
into a turbulent flow, they quickly distribute themselves 
throughout the volume of the fluid; further stirring leaves 
the particle density distribution uniform. The bouyant par- 
ticles used in this experiment have a density much less than 
that of the fluid on which they move. Hence their behav- 
ior is entirely different. Now, their motion is confined to 
the surface of the turbulent fluid. If the spatial distribu- 
tion is initially uniform, at say, t = 0, continuous stirring 
of the underlying fluid evolves this spatial distribution into 
string-like structures. Ultimately a steady state is reached, 
at which time the floaters occupy a fractal dimension much 
less than 2 [14] . This time evolution into a fractal is a generic 
effect. It occurs even if the interactions between the floaters 
is negligibly small. A common manifestation of this phe- 
nomenon is the coagulation of scum on the surface of the 



sea, as is often seen in an ocean harbor. 

The transient evolution of the floaters is studied by uni- 
formly distributing the floaters at t = (as discussed in 
2.). At times t > the particles are subject to the underly- 
ing turbulent flow until their spatial distribution reaches a 
steady state. To study the floater's fractal distribution, the 
spectrum of fractal dimensions D q (t) is measured at subse- 
quent times t during the transient evolution. One may ar- 
gue that the generalized fractal dimensions D q of any sys- 
tem are meaningful only in the limit of the evolution time 
t — > oo. However, there exist instances of fractal objects re- 
alizable in real or configuration space (such as shear flows 
[35], cement gels [31], and [17]) that are amenable to inves- 
tigation of their transient behavior. 

It is essential to realize that the floaters in this exper- 
ment are passively advected by the underlying flow. They 
are small enough to follow the velocity field of the turbulent 
sea on which they move, in a plane that has coordinates x, 
y, z = 0. There are, of course, waves on the surface, and 
they can drive the motion of the floaters. Separate studies 
have shown that the amplitude of the surface waves is small 
enough to have a negligible effect on the particle motion 
[18]. Though the water molecules can have a downward ve- 
locity component at all values of the depth z, the floaters 
cannot follow them. That is why they coagulate and dis- 
perse in the plane of their motion x, y with z = 0. In this 
sense the floaters form a strongly compressible system. Us- 
ing the following definition of the dimensionless compress- 
ibility, 
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With this definition, C must lie between zero (incompress- 
ible) and 1 (potential) for an isotropic flow field. Experi- 
mentally, C is close to 0.5 [18]. 

The coagulation phenomenon described above was 
demonstrated and analyzed by J. Sommerer and E. Ott 
(S&O) [16]. The solution was rather gently stirred via a 
pulsing jet, its motion being slow enough that they could 
measure the steady-state fractal dimensions of the surface 
particles and the two Lyapunov exponents as well. The lat- 
ter parameters, Ai > A2 define the rate at which initially 
close particle pairs separate in time. Since the pattern be- 
comes string-like in the steady state, the largest exponent 
Ai is positive, and the other one A2 is negative. Since the 
total area ultimately occupied by the floaters decreases, 
Ai + A2 < 0. The information dimension (Z?i) of the frac- 
tal pattern can be related to the dynamics of the system 
through the Kaplan- Yorke dimension Dky = 1 + Ai/ 1 A2 1 
[21]. For a two dimensional system, Dky=D\. In a sepa- 
rate study [17] (but with a similar experimental setup as in 
[16]), the correlation dimension D2 of the floating particles 
was measured, both at an initial time t ~ when D2 — 2 
and in the steady-state where 1 < D2 < 2. In that exper- 
iment, the correlation dimension could not be measured 
during the transient state because of poor scaling [17]. In 
the present experiment, we observe a robust scaling of the 



correlation sum 62(7") during the transient evolution, see 
Figure 3, from an initial state 1)2(0) — 2 to a steady-state 
value 1)2(00) ~ 1.25. 

The present experiment differs from that of S&O in that 
the stirred fluid, water, could be driven into a strongly tur- 
bulent state. The Taylor microscale Reynolds number, Re\ 
is approximately 160 (see Table I). This large Re\ estab- 
lishes a well defined inertial range of the flow (sec [18]). 
For a more thorough exploration of this type of flow, see 
[18,20]. Since the experiment of S&O, there have been sev- 
eral theoretical advances pertaining to clustering phenom- 
ena in turbulent flows. These theories utilize the statistics 
of stretching rates below the dissipative scale 77 of turbu- 
lence (Table I) to predict a multi-fractal particle distribu- 
tion for compressible flows [26]. The measurements made 
in this experiment are strictly for scales greater than the 
dissipative scale 77 (see section 2), where there is currently 
no theoretical guidance [26,22]. 

Before describing the experiment in detail, it is helpful 
to observe the pattern of the floaters at a sufficiently long 
time so that the steady state has almost been achieved. 
Figure 1 shows the distribution of the Lagrangian particles 
(discussed below) at t = 0.15 s and t = 1.5 s, their spatial 
distribution being uniform at t = 0. Were an image made 
at t =0, the particle distribution (shown as dots) would be 
uniform. The blank white spaces are due to the finite initial 
seeding procedure and do not affect the results. 
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Fig. 1. Visualization of particles clustering from a uniform initial 
distribution of 10 5 tracer particles for Re,n=169 at t=0.15s for (top) 
and t=1.5s (bottom). By t=1.5s the particle distribution has reached 
a steady-state, which occurs within several Large Eddy Turnover 
Times r (Table 1.) 



2. Experiment 

The 1 m x 1 m tank is filled with water to a height of 
30cm. The tank is large compared to the camera's field of 
view. The turbulence is generated by a large pump con- 
nected to a network of rotating jets in a plane 10 cm above 
the tank floor. See Figure 2 for a schematic of the experi- 
mental setup. The arrangement creates uniform turbulence 
in the center of the tank, and also moves the source of tur- 
bulent injection far from the fluid surface where the mea- 
surements are made [18]. With this scheme, surface waves, 
which cannot be avoided, do not exceed an amplitude of 
~1 mm [18]. It is necessary that the surface of the tank 
be freshly cleaned before each set of measurements. Oth- 
erwise, amphiphiles form a continuous layer on the surface 
and prevents the floaters from moving freely under the ac- 
tion of the turbulence [18]. 

The hydrophillic particles chosen here are subject to cap- 
illary forces which are very small compared to forces com- 
ing from the turbulence, and do not affect the results as 
they do in [36,37]. The non-inertial character of the parti- 
cles is minimal because the Stokes number St is small: St = 
T s v rm s/ & — 0.01, where a is the particle radius, v rms is the 
RMS velocity of the turbulent fluid at the free-surface, and 
t s is the stopping time of the particle [38] . 

During an experimental run, the floating particles (50/xra 
diameter and specific gravity of 0.25) are constantly seeded 
into the fluid from the tank floor, where they undergo tur- 
bulent mixing as they rise due to buoyancy and are uni- 
formly dispersed by the time they rise to the surface. Once 
at the free-surface, their motion is constrained to the two- 
dimensional surface plane. Their motion is tracked via a 
high-speed camera (Phantom v. 5) situated above the tank. 
The camera field-of-view is a square area of side length L = 
9cm. The constant particle injection is necessary to replace 
particles at the surface during the experiment. The source 
and sink structure at the surface fluctuates in both time 
and space, which can cause particles to leave the camera's 
field of view. 

Instantaneous velocity fields are measured using an in- 
house developed particle imaging velocimetry (PIV) pro- 
gram which processes the recorded images of the floaters. 
The constant injection of particles ensures that surface 
sources and sinks receive an adequate coverage of particles 
on the surface. The local particle density at the surface de- 
termines the average spacing of the velocity vector fields 
produced by the PIV program. The resulting velocity vec- 
tors are spaced (on average) by Sx = 2.5 rj over both sources 
and sinks. This vector grid spacing is important for the La- 
grangian particle evolution scheme, which is discussed be- 
low. The camera's height above the water surface was cho- 
sen so that a pixel size is roughly 0.1mm, comparable to 
the dissipative scale of the turbulence. 

The measured velocity field was then used to solve the 
equation of motion for Lagrangian particles : 



(2) 



where v(xj, t) is the velocity field and Xi = (ar», are the 
individual particle positions. To achieve accurate results 
for the Lagrangian particle evolution, the vector fields used 
in Eq. 2 were interpolated from the experimentally deter- 
mined velocity vectors via a bi-cubic interpolation scheme 
developed for numerical simulations, as discussed in [27] 
and implemented in [18]. This scheme uses the smooth flow 
between grid points separated by length scales compara- 
ble to r\ to interpolate the velocity field between measured 
velocity grid points. To use this scheme it is necessary for 
the measured velocity grid spacing to satisfy the criterion 
Sx < 7r?7, where 5x is the above mentioned average mea- 
sured velocity grid spacing. We have tested to ensure that 
the results do not depend on the velocity grid spacing by 
varying the spacing from 8x — 2.5 77 to Sx = 4 rj. 

The Lagrangian particle tracks evolved by Eq. (2) are 
then used for the measurements presented in this work. 
This is the method used to achieve a uniform distribution 
of floaters at t = 0. Visualization of these Lagrangian trac- 
ers can be seen in Figure 1. The experimental setup is dis- 
cussed in more detail in [18]. Data were taken for several 
values of Re\ ~ 150 — 170 with an average Re\ ~ 160. 
Since the measurements show no systematic variation with 
the Re\ over this range, each experimental run was aver- 
aged to decrease measurement errors. Turbulent parame- 
ters measured at the surface are listed in Table 1. All of the 
statistics presented below were obtained by evolving ~ 10 5 
Lagrangian particles per frame. Tests were performed to 
ensure that the number of tracers provided adequate statis- 
tics to calculate D q (t). The initial homogeneous seeding of 
particles was varied from 10 5 to 4 x 10 5 and the results were 
insensitive to this variation. Since computation time goes 
roughly as N 2 because of Eq. (4), no more than 10 5 tracers 
were used. 

3. Results and Analysis 

We investigate the inhomogeneous particle distributions 
by measuring their time-evolving fractal dimensions. Frac- 
tal dimensions are mathematical representations of com- 
plex patterns and provide measures of spatial (or temporal) 
dependence at a variety of scales. For any q, the spectrum 
of fractal dimensions is [19]: 
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where the correlation functions C q {r) are defined: 
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For q = 2, the quantity contained inside the brackets in 
Eq. (4) is the probability of two randomly chosen points 
(here our passive tracers) being within a distance r of one 
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Fig. 2. Schematic of the top-view (top panel) and side-view (bottom 
panel) of the experimental setup. 36 rotating capped jets are placed 
horizontally on the tank floor (shown as randomly oriented Z-shaped 
patterns) that pump water into the tank re-circulated by a 8hp 
pump. The central region of the water surface is illuminated by a 
laser-sheet. A high-speed digital camera suspended vertically above 
this central region captures images of the light scattered by buoyant 
particles (50 fim hollow-glass spheres of specific gravity 0.25). 

another. For q > 3, the bracketed function is the number 
of q-tuplets of points (particles) whose pairwise distance is 
less than r. Here N are the total number of tracer particles 
(forced to be constant), 6 is the heaviside step function, 
and Tij is the distance between particle i and j. This algo- 
rithm for determining the spectrum of fractal dimensions 
was given by Hentschel and Proccacia [15]. 



Table 1 

Turbulent parameters measured at the surface. Measurements are 
made at several values of the Re\ with an average Re\ ~ 160. The 
parameters listed are averages, with deviations less than 10%. 



Parameter 


Symbol used in text 


Measured 
Value 


Taylor microscale A (cm) 


A \l ((S^/S^) 2 ) 


0.47 


Taylor Re^ 


R v rms \ 


160 


Integral Scale lo (cm) 


, f dr hi(^)°iif)) 
J <(«„(«»') 


1.42 


Large Eddy Turnover time 
(LETT) T (s) 




0.43 


Dissipation Rate £diss 
(cm2/s3) 

Kolmogorov scale rj (cm) 


^ ss =io,((^) 2 ) 


6.05 
0.02 


RMS Velocity vrms (cm/s) 


Vrms = \/ (v 2 ) ~ (v) 2 


3.3 


Compressibility C 


-<(;,,)*)■ 


0.49 ± 2% 



To calculate D q , the log of the correlation sum (Eq. (4)) 
is plotted versus the log of r. The range of r over which the 
plot is a straight line is the scale-free (or scaling) region. 
The slope of the line dlog (C g (r ))/dlog(r) is the value of D q . 
The abscissa is r in units of the dissipative scale r\ (Table 1). 

At t — the particles are uniformly distributed, D q (Q) = 
2. Figure 3 is a plot of log(C2(r)) versus log(r) for t = 
and subsequent times. The scaling range is seen to span 
the interval 10 5 < r/rj < 10 2 . The ratio of the integral 
scale lo to the dissipative scale ?y is lo/r/ ~ 70. The initial 
homogeneous particle separation is roughly 7/. At t = Q, 
D-2.it = 0) ~ 2 for 3 < r/77 < 70. For scales r/rj < 2, D 2 (t = 
0) ~ 0, indicating the point-like nature of the Lagrangian 
tracer particles at scales less than initial tracer spacings. 
Thus, to ensure that the D q are defined such that D q > 
for all times t, the only measurements presented here are 
for r/rj > 2. 

Figure 4 shows measurements of D q for Re\ ~ 160 for a 
range of q = 1 to 6 at various times. The measurements at 
t = 1.5s are in the steady state (discussed below). The re- 
sults indicate a multi-fractal distribution over what is usu- 
ally considered to be the inertial range of the flow. How- 
ever, other experiments [39,40] observe that the transition 
from the inertial to dissipative range occurs at spatial scales 
greater than 77. Because of this, it is difficult to label the 
measurements here as being strictly inertial or dissipative. 
While the theory in [26] predicts a multi- fractal particle dis- 
tribution for compressible flows below the dissipative scale, 
no theory exists for the inertial range. The results in this 
work do agree qualitatively with the existence of inertial- 



4 



range scaling seen in a numerical simulation [30]. The fol- 
lowing analysis of the time-evolution of the dimensions D q 
will focus on the so-called correlation dimension D2 ■ 




Fig. 3. Correlation functions C%{r) averaged over several Re\ ~ 160 
for various times in the experiment. The interval of r that exhibits 
scale-free behavior is r/rj ~ 10 2 to 10 ' 5 . The inset shows the value 
Z?2 calculated from the main figure. Note that these scales are in the 
inertial range of the flow (Table I). 
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Fig. 4. Range of values D q measured for Re\ = 169 for various times 
in the experiment. These results indicate a multi-fractal distribution 
of the particles for the inertial range of the flow after the particle 
distribution evolves from a homogeneous distribution at t = s. 

The experiment was performed at several Rex — 150 — 
170, with an average Re\ ~ 160. From each of these exper- 
iments, 17 time traces of the various D q (t) were measured 
starting with an initial homogeneous distribution. These 
individual time traces were then used to produce an ensem- 
ble measurement of D q (i), which is analyzed subsequently. 
The total time spanned by the ensemble of experiments is 
approximately 25s, or roughly 60 LETT's. For t > 1.5, all 
of the D q measured fluctuate around a steady-state limiting 
value D q (oo). For D2(t), 1)2(00) = 1-25. Figure 4 shows an 
example of D q (t) for Re\ — 169 measured at various times 
t < 1.5s, after which the dimensions saturate to values of 
D q {oo). These values D q (oo) are approximately those at 
D q {t = 1.5s). The error bars reported in Figure 4 are due 
to the error in the fits used to obtain the D q (t). Similar 
errors are found from the statistical ensemble averaging of 
the D q (t). 



Our main result appears in Figure 5. This figure shows 
that D2(t) (ensemble averaged) decays at an exponential 
rate from 1)2(0) = 2 to 1)2(00) : D2(t)— 1)2(00) oc exp(-jt), 
where 7 = 2.4 ± 0.1Hz. The decay time r = I/7 is 0.4s, 
which is approximately one large-eddy turnover time To 
(Table 1 ) . This is the typical time for the largest " eddies" to 
significantly distort in a turbulent flow [24]. Several other 
D q (t) also experience an exponential decay from an initially 
homogeneous state D q (0) = 2 to their steady-state limiting 
value Dq(oo). The inset of Figure 5 shows the exponent 7(g) 
characterizing the exponential decay. One can see that 7(9) 
is approximately independent of q for q — 1..3. However, 
for q > 3 the error in the measurement makes it difficult to 
predict that for large q, D q (t) will decay exponentially. It 
may also be worth noting that the value of 7(17 = 1..3) is 
close to the value of the smallest Lyapunov exponent (but 
the largest in magnitude) measured in [23] (A2 ~ —2Hz). 
It remains to understand why this decay has exponential 
form. 




t(s) 

Fig. 5. D2(t) averaged over several Re\ ~ 160. The result is an ex- 
ponential decay of D2(t) from D2(0) = 2 to roughly 1.25 in approxi- 
mately 1.5 s. The solid line is a best fit to the data. The inset shows 
the decay exponent 7(g) (discussed in text) as it varies with q. 

4. Summary 

A compressible system of particles in free-surface turbu- 
lence represents an instance of a chaotic attractor in real 
space. We study the time evolution of the low-order frac- 
tal dimensions D q for an ensemble of particles floating on 
a turbulent tank of water with the initial value of D q = 2 
(17=1. ..6). The system evolves in a time of the order of the 
lifetime of the largest eddies, to a steady state where the 
measured D q (t) approach a value that is slightly greater 
than 1, implying the formation of string-like structures. The 
correlation dimension D 2 (t) evolves exponentially as the 
steady-state is approached. It is not possible, so far, to de- 
duce these observations from the Navier-Stokes equations 
in the inertial range of the compressible flow studied here. 
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